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>• 1 Abstract 

, We explore the implications of possible neutrino oscillations, as in- 

dicated by the solar and atmospheric neutrino experiments, for the 
■ cold plus hot dark matter scenario of large scale structure formation. 

I/"") , We find that there are essentially three distinct schemes that can ac- 

commodate the oscillation data and which also allow for dark matter 
neutrinos. These include (i) three nearly degenerate (in mass) neu- 
trinos, (ii) non-degenerate masses with v r in the eV range, and (iii) 
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nearly degenerate — v T pair (in the eV range), with the additional 
■ possibility that the electron neutrino is cosmologically significant. The 

Ci ■ last two schemes invoke a 'sterile' neutrino which is light (<> eV). We 

discuss the implications of these schemes for — v e and — v T oscil- 
lation, and find that scheme (ii) in particular, predicts them to be in 
the observable range. As far as structure formation is concerned, we 
compare the one neutrino flavor case with a variety of other possibil- 
ities, including two and three degenerate neutrino flavors. We show, 
both analytically and numerically, the effects of these neutrino mass 
scenarios on the amplitude of cosmological density fluctuations. With 
a Hubble constant of 50 km s _1 Mpc -1 , a spectral index of unity, and 
^baryon = 0.05, the two and three flavor scenarios fit the observational 
data marginally better than the single flavor scheme. However, taking 
account of the uncertainties in these parameters, we show that it is 
premature to pick a clear winner. 
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1 Introduction 



The idea that the dark matter in the universe may contain both a 'cold' as 
well as a 'hot' component M first arose as a serious possibility within the 
framework of grand unified theories (which also inspired inflation, cosmic 
strings, baryogenesis, etc.). In the original models the cold component was 
the axion, although nowadays the lightest supersymmetric particle (LSP) 
often plays that role. Lightly massive neutrinos are the hot component, 
and it was usually assumed that the neutrinos would display the hierarchical 
mass spectrum characteristic of the quarks and charged leptons, with a single 
neutrino (e.g., u T ) dominating the mass density of the hot component. 

However, it was noted from the beginning jl| that because of their special 
nature, nearly degenerate massive neutrinos is a logical, albeit not necessar- 
ily the simplest, possibility. The cosmological implications of neutrinos that 
are closely degenerate in mass has received much recent attention || f|, |], 
inspired to some extent by the remarkable series of neutrino oscillation ex- 
periments (solar, atmospheric, and more recently accelerator) which suggest 
that two or more neutrino flavors may contribute to the hot dark matter in 
the universe. 

Before discussing the impact of neutrino oscillation experiments on the 
cold plus hot dark matter (C+HDM) scenario, let us briefly recall why the 
latter has attracted so much recent attention. From the mid to the late 
1980s evidence was mounting that the so-called "standard" cold dark matter 
(CDM) scenario had trouble providing a consistent explanation of small (~ 
galactic) as well as large scale structure. On the contrary, a C+HDM model 
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(with Q u ~ 0.15 — 0.35) provided a far more consistent fit to the data ||. In 
order to provide additional tests, the quadrupole anisotropy of the microwave 
background expected in C+HDM models was estimated in 1989 Q, |8| and 
compared with the CDM prediction. Normalized to the 'small' scale data, it 
was found that the C+HDM prediction exceeded the CDM value by a factor 
of about two, and this was dramatically verified when the COBE team II 
made its announcement in 1992 ||10| . 

In the last two years, particularly after extensive numerical simulations of 



galaxy and cluster formation JTTj] and more detailed analytical work [12, [13 



additional support for the C+HDM model has emerged. It appears that this 
model provides the simplest (consistent) realization of an inflationary sce- 
nario for large scale structure formation. It has recently been pointed out 
[||. f|], (and our results are in agreement) that with two or more neutrino fla- 
vors contributing to the hot dark matter, somewhat better fits to the present 
data may be possible. For fixed VL U there is also some dependence both on the 
spectral index as well as the Hubble constant. It must be admitted, though, 
that within the framework of grand unified theories, it may not be easy to 
realize scenarios with two or three nearly degenerate (in mass) neutrinos. 
Furthermore, in some of the scenarios we consider, one needs to invoke a 
fourth (sterile) neutrino which is nearly degenerate in mass with the electron 
neutrino. Precisely why this singlet state can be so light is an important 
question for future research. 

The layout of the paper is as follows. In section 2.1 we list three sets of 
observations related to solar, atmospheric, and dark matter neutrinos which 
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provide hints for non-zero neutrino masses. For completeness, we also sum- 
marize the results of the LSND experiment at Los Alamos, although we will 
not include them among the constraints to be satisfied. We then proceed 
to a discussion of three distinct schemes which are consistent with the ob- 
served solar and atmospheric v oscillations, and which also allow for neutrino 
dark matter. We discuss the implications of these schemes for a variety of 
experiments, particularly the P M — v e and — v r oscillation experiments at 
accelerators. In section 3 we address the main issues raised by neutrino oscil- 
lations for cosmological structure formation within the C+HDM framework. 
We review the physics of the differences between having one and more than 
one massive neutrino flavor, giving some analytic formula to estimate the 
scale (section 3.1) and amplitude (section 3.2) of the effects on the growth of 
primordial density fluctuations. In 3.3 we give results of some more detailed 
numerical computations of the cosmological models based on the solutions to 
the oscillation data, including some models, in 3.4, motivated by the LSND 
results. Finally, in section 3.5 we discuss the dependence on the cosmolog- 
ical parameters (Hubble constant, spectral index n, and baryonic density 
fraction) of conclusions drawn from comparing specific models to large scale 
structure data. 

2 Neutrino Masses and Mixings 
2.1 Observational Hints 

One of the greatest challenges in particle physics today pertains to the issue 
of the neutrino masses. There are indications from a variety of experiments 
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that one (or more) of the known neutrinos may possess non-zero mass(es). 
In this section we briefly summarize the relevant observations based on solar 
and atmospheric neutrinos, as well as present arguments for neutrino dark 
matter. [For completeness, we also summarize the recent findings of the 
LSND experiment at Los Alamos.] In section 2.2, we discuss scenarios for 
neutrino masses and mixings which can simultaneously accommodate these 
observations. We find that there are essentially three viable schemes, and 
remarkably enough, all three will be tested in ongoing and planned neutrino 
oscillation experiments (particularly Los Alamos LSND, Rutherford KAR- 
MEN, the CERN CHORUS/NOMAD and Fermilab E803 experiments). 

1. Solar neutrino puzzle: The apparent deficit in the flux of solar z/ e 's 



14 1 that has persisted in the Chlorine experiment for over two decades 



15 1 has been confirmed in the last several years by three independent 



experiments: the Kamiokande water Cerenkov detector and the 



SAGE |HJ and GALLEX [|THJ radiochemical experiments. These four 
observations, which probe different energy regimes in the solar neutrino 
spectrum, can all be simultaneously explained in terms of two flavor 
neutrino oscillations. The matter enhanced MSW |19| oscillation for a 
two flavor {y e — v x ) system admits two branches [pOfl : (i) The small angle 
non-adiabatic solution, which requires sin 2 26 l = (3.5 x 10~ 3 to 1.5 x 
1(T 2 ) and Am 2 = (m 2 -m 2 ) = (3.4 x 1(T 6 to 1.2 x 1(T 5 ) eV 2 . (ii) The 
large angle MSW solution requiring sin 2 26' = (0.6 — 0.9) and Am 2 = 
(7 x 10~ 6 to 5 x 10 -5 ) eV 2 . Here v x could be u^, v T or a sterile neutrino 
v s for the small angle MSW solution, while v x = (z/ M or v T ) in the large 



angle MSW solution. There is also the {y e — v x ) vacuum oscillation 



solution [2l] for v x = (y^ or u T ), with sin 26 = (0.6 — 1.0) and Am ~ 
10- 10 eV 2 . 

2. Atmospheric neutrino anomaly: The ratio of muon to electron neu- 
trinos from atmospheric cascades measured by the Kamiokande 
and the 1MB |£| experiments appear to suggest a deficit of about 
a factor of 2 when compared to Monte Carlo simulations || . For 



sub-GeV neutrinos (i.e., neutrinos with energy less than ~ 1 GeV), 
Kamiokande finds for this ratio of ratios (which is expected to be 1 in 
the absence of neutrino oscillations) R = 0.60^'^ ± 0.05 which is in 
good agreement with the 1MB value of R = 0.54 ±0.05 ±0.12. Recently 
Kamiokande has also reported results for the multi-GeV neutrinos p5f 
where the ratio is R = 0.57±g;g| ± 0.07, in agreement with the sub- 
GeV data. This apparent deficit of u^s or the excess of z/ e 's can be 
attributed to either [y^ — v^) oscillations with sin 2 26' = (0.65 — 1.0) and 
Am 2 = (5 x 10 -3 to 3 x 10~ 2 ) eV 2 , or to [y e — v^) oscillations with 
sin 2 2# = (0.55-1.0) and Am 2 = (7x 10~ 3 to 7x HT 2 ) eV 2 . The resolu- 
tion of the anomaly in terms of (v^ — v s ) oscillations, where v s is a sterile 
neutrino, will run into difficulty with primordial nucleosynthesis calcu- 
lations which require that |2£| Am 2 sin 2 2# < 1.6 x 10~ 6 eV 2 . (Note that 
the solar (v e — v s ) oscillation parameters satisfy this constraint.) Among 
the other atmospheric neutrino experiments, only Soudan II |27j sees an 
anomaly with a preliminary value oi R = 0.69±0.19±0.09. The Frejus 
p8| results are marginally in conflict with the combined Kamiokande 
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and 1MB data. Considering the low statistics of Frejus experiment 
compared with the high statistical significance of Kamiokande and 1MB 
data, we shall ignore this marginal discrepancy in our theoretical dis- 
cussions. 

3. Hot component of dark matter: As already noted in the introduction, 
and further explained in the next section, a combined fit to the COBE 
data as well as the data on distribution of galaxies on the large and 
small angular scales is difficult to achieve within the cold dark matter 
scenario. The simplest consistent scenario requires a significant hot 



component (15-30)% in the dark matter [y| |29|, leading to the 
C+HDM scheme, with neutrinos being the natural candidates for the 
hot component. The mass of the neutrino comprising the hot dark 
matter should be in the few eV range. 

4. 17^ — V e oscillation at Accelerators: Recently the LSND experiment at 
the LAMPF facility in Los Alamos has reported positive evidence for 
— V e oscillations |JT] . If these results survive further scrutiny and 
there is independent verification (e.g, at the KARMEN experiment 
fl32|), it certainly will have a strong impact on particle and nuclear 
physics as well as cosmology and astrophysics. The initial LSND data, 
if interpreted in terms of two-neutrino oscillations suggest a mixing 
probability of Pu^v e = (0.341^^0.07)%, when oscillation constraints 
from KARMEN and BNL-E776 |33] accelerator searches as well as re- 



actor neutrino constraints from Bugey facility |34| are folded in. In 
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the discussion that follows, we shall take into considerations the sen- 
sitivity of the LSND and other accelerator and reactor experiments 
in searching for 77 M — V e oscillations, but we shall not demand that 
the theoretical scenarios discussed correctly reproduce the precise nu- 
merical values reported by the LSND collaboration. We feel that this 
position is justified at this time, especially in view of the fact that a 
different analysis of the LSND data has yielded considerably weaker 



limits on the mixing parameters [35|. We find it remarkable however, 



that independent of the LSND data, two of the three scenarios which 
can accommodate points (l)-(3) above, imply an observable signal at 
accelerator (y^ — V e ) oscillation experiments. 

2.2 Theoretical Schemes: 

Let us focus on the solar, atmospheric and hot dark matter neutrinos, 
leaving out the LSND results which will not play a role in determining the 
viable scenarios. It is clear that the mass-splittings required for explaining 
these three observations do not overlap (Am 2 ~ 1CT 5 eV 2 for solar neutrinos, 
Am 2 ~ 10~ 2 eV 2 for atmospheric neutrinos and m u ~ (1 to few) eV for hot 
dark matter). Assuming that there are no additional light neutrinos, one 
concludes that the three neutrinos (u e , u T ) should be nearly degenerate in 
mass. (Note that we often identify a mass eigenstate through its dominant 
flavor. The respective masses are mi, m^, and 777.3.) This leads us to our 
first scenario. 

Scenario (i): Three nearly degenerate neutrinos: Assume that v e , and v T 
have a nearly common mass of about (1 to 3) eV. Their masses are split 



by small amounts, such that Am^ = (m| — mf) ~ 10~ 5 eV 2 and |Am| 3 | = 
|(m| — ~ 10~ 2 eF 2 . Phenomenological neutrino mass matrices that are 
consistent with these assumptions are easily constructed |51j|| . This scheme 
will account for the solar neutrino data via iv e — v ^ ) MSW oscillations and the 
atmospheric neutrino anomaly in terms of (u/j, — u T ) oscillations. The required 
mixing angles are free parameters and can be adjusted to the suggested 
values. The {y e — v T ) mixing angle can be arbitrarily small, but it can also 
be as large as about 0.05. According to this scenario, neutrinoless double 
beta decay (PPou) is at (or even above?) the present experimental limit 
(modulo nuclear matrix element uncertainties), and so that it can be ruled 
out if PPq u is not observed in the near future. 

Also according to this scheme the — v T oscillation experiments at ac- 
celerators (CHORUS/NOMAD, E803 |37|, |38|), as well as v e - oscillation 
searches at reactors, should fail to find anything interesting. In both cases 
the relevant Am 2 < 10~ 2 eV 2 , which is not in the range accessible in these 
experiments. However, this statement is strictly true only to the extent that 
the leptonic mixing matrix is assumed to be unitary. Deviation from uni- 
tarity can produce oscillation signals, although there should be no spatial 
dependence. This situation can arise if the ordinary neutrinos mix with a 
neutral heavy lepton with mass on the order of the electroweak scale. Let us 
consider this in more detail. 

For definiteness, let us assume that the three light neutrinos have an 
admixture of a heavy fourth generation neutrino with mass greater than 
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M z . Writing 

4 

Ve = ^ZU ei Vi (1) 
i=l 

and similarly for the weak eigenstates u^, u T , one sees that muon neutrinos 
produced in 7r decays will be the weak eigenstates with the heavy z/ 4 compo- 
nent removed. The v e and states are no longer orthogonal, leading to an 
apparent oscillation probability given by [^] 

= i-o^mi 2 • (2) 

Here we have set the Am 2 s to zero since they are all < 10~ 2 eV 2 and thus 
negligible for accelerator neutrino experiments. If the mixing parameters 
U e 4, ~ (j to |), then the oscillation probability at the LSND experiment 
will be in the range of (2 — 3) x 10~ 3 with no position dependence. The 
data presented in Ref. j^TJ is consistent with such an interpretation, as 
can be inferred from the large Am 2 region of their Fig. 3. Various other 
constraints on such leakage to a heavy fourth generation neutrino has been 



studied in |40|, where it is shown that the best existing limit is from neutrino 
oscillation experiments, — v T oscillation probability of (2 — 3) x 10~ 3 can 
also be achieved by a similar mechanism, but this may be more difficult to 
measure. 

If the heavy neutrino is a standard model singlet, there are stringent limits 
from Z-boson decays which make the non-unitary oscillations unobservable 
at accelerators. To see this, consider the mixing of ordinary neutrinos with 
a heavy isospin singlet neutral lepton N of mass m. In the presence of such 
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mixings, the invisible decay width of Z will be modified so that the number 
of effective neutrino species coupling to Z is given by 



AT - 3+ (1- \a 



2^ 



-(l + H 2 )+2F|a| 2 + F'(l- \a\ 2 ) (3) 



where F and F' are phase space factors given by 

F = 1 - -x 2 + -x 3 
2 2 



F' = Vl -4x 2 (l -x 2 ) (4) 

with x = m/mz- It is to be understood here that for m < mz/2, all terms 
in eq. (3) will contribute to N u , while for m > mz, the F and i 7 " terms 
do not contribute. If mz/2 < m < mz, the term proportional to F' should 
be set to zero, while F is non-zero. The mixing parameter \a\ is equal to 
the (4,4) element of the unitary matrix U v which diagonalizes the neutrino 
mass matrix. It is related to the leptonic mixing matrix U appearing in the 
charged current via 

H 2 = i-\u ei \ 2 -\u^\ 2 -\u Ti \ 2 . (5) 

Although U, now a 3 x 4 matrix is not unitary, the u e — transition 
probability is still given by eq. (2), since each row of U is normalized to 
1. It becomes apparent from eq. (2) and (5) that P v ^- Vli < (1 — |a| 2 ) 2 /4. 
If the mass of A^ is greater than mz, then F = F' = 0, so that P Vf ,^ v ~ 
[(3-A" ; ,)/(l + |a| 2 )] 2 /4 = [AA"^/(l + |a| 2 )] 2 /4. AN V can be as large as 0.042 at 



1 sigma I^TJ, which when combined with | ct | 2 ~ 1 yields P Ue -u — 1-1 x 10~ 4 , 
too small to be observable. (Even with AA"^ = 0.067, the 2 sigma value, 
Pv e -u — 2.8 x 10 -4 .) If the mass of A^ is less than mz/2, then both the F and 
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F' terms in eq. (3) will be relevant, making N v more consistent with 3. (As 
F,F' — > 1, N v — > 3.) However, such a scenario is ruled out by direct search 
limits for a heavy neutral lepton decaying into the usual charged leptons 
(e, fi), as should happen here. (The limit on the mass of such a neutral 
lepton is m > 46 GeV).) If the mass of iV obeys m z /2 < m < m Z) then 
F^O, but F' = in eq. (3). In this case, noting that F is at most 11/16 (for 
x = 1/2), we see that Pve-v^ can be as large as (4AN u /5) 2 ~ 1CT 3 , which 
is close to present experimental sensitivity. However, since the branching 
ratio Br[Z — > (Nv + Nv)] ~ 2 x 1CT 4 , few hundred such events should have 
been observed at LEP with the N decaying subsequently into This 
possibility is also excluded based on the non-observation of such events. 

We conclude that the case of isosinglet neutrino does not lead to any 
observable deviation from unitarity in the accelerator neutrino oscillation 
experiments. Note that there is no such constraint on a sequential fourth 
generation neutrino (or neutrinos in vector-like families). 
Scenario (ii): No degenerate neutrinos: If there is no degeneracy in the neu- 
trino masses and we assume some kind of a mass hierarchy, accounting for 
the solar, atmospheric and hot dark matter neutrinos would require the in- 
troduction of a light (sterile) neutrino v s . There is then just one consistent 
scheme in this case. As far as we are aware, this possibility has not been 
discussed in the literature, so we shall elaborate on it. 

By assumption, in this scenario the known neutrinos are non-degenerate, 
and indeed have hierarchical masses. The solar neutrino puzzle is resolved via 
(v e — v s ) MSW oscillations such that m 1 Vs —m 2 Ue ~ 1CT 5 eV 2 . The atmospheric 
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neutrino anomaly is explained via (y^ — v e ) oscillations with m v ~ 10 _1 eV. 
The tau neutrino with a few eV mass constitutes the 'hot' component of dark 
matter. 

As far as (z/ M — v e ) oscillations at accelerators are concerned, the "direct" 
transition is not possible because Am 2 ~ 10~ 2 eV 2 . However, it has recently 
been pointed out EI that "indirect" (y^ — V e ) transition via a virtual V T in 
the few eV mass range can still occur in such a scheme. For this to be in 
the experimental reach, both [y e — u T ) oscillations (at reactors) and (y^ — u T ) 
oscillations (at accelerators) should be in the observable range. As shown in 



[ fi^j , the present limits from Bugey reactor [33| on u e — v T oscillations, and 



the CHARM-II, Fermilab E531 and CDHS limits [|| on (p^-p t ) oscillations 
translates into an observable "indirect" (z/ e — u^) oscillations for a v T mass 
in the few eV range, which can account for the LSND results. 

Let us point out that in this second scenario, the allowed MSW parameter 
region will be somewhat shifted because of the rapid (y e — v^) oscillations 
from the sun to earth, leading to a further depletion of u^s by a factor of 
|sin 2 2^ eM ~ (| to \). This prediction will be tested in the forthcoming solar 
neutrino experiments. 

This scheme may be on a somewhat better theoretical footing than the 
first one because no mass degeneracy is assumed. A light 'sterile' neutrino 



state, following [44], may arise from the hidden sector of some fundamental 
theory. If the hidden sector contains SU(3) x SU(2) x £7(1), one can define 
an unbroken 'parity' at the Lagrangian level. The vrs needed for the see- 
saw mechanism, having no gauge quantum numbers can freely mix from the 
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observable and hidden sectors. This would lead to mixing between the light 
neutrinos and their mirror partners (which are also light). Gravitationally 
induced interactions could also mix the neutrinos from the two sectors. 
Scenario (Hi): Nearly degenerate [y^ — u T ) pair: This scheme also requires a 
light sterile neutrino v s . In addition, an approximate lepton number symme- 
try such as L e + — L T is necessary to make and v T nearly degenerate in 



mass |36| , |44] , |45| . When this symmetry breaks by a small amount, maximal 
(z/ M — v T ) mixing occurs, facilitating a resolution of the atmospheric neutrino 
problem. The solar neutrino puzzle is explained via iy e — u s ) oscillations. 
The common mass of and v T is assumed to be ~ few eV, so that they 
constitute a 'two-flavor' hot dark matter. For accelerator neutrinos, there is 
direct [y e — v^) oscillations at an observable level. In particular, the LSND 
data can be accommodated. However, [v^ — u T ) oscillations will be beyond 
the reach of CHORUS /NOMAD and E803. 

Aside from the two-flavor hot dark matter component, this scenario has 
other special cases relevant for cosmology. The solar neutrino data requires 
the mass splitting Am 2 between v e and v s to be ~ 10 -5 eV 2 , but the masses 
themselves may be in the cosmologically interesting range of a few eV. Since 
this pair behaves like Dirac neutrinos, there is no conflict with /3/3o v con- 
straints. For cosmology however, the hot dark matter may be a linear com- 
bination of v e and [v^ + u T ). If the total mass of (z/ M , u T ) pair exceeds that 
of v e , we call it a (2+1) C+HDM scheme, and if the v e mass dominates over 
(v^, u T ), it is a (1+2) scheme. Of course, one could also recover the one 
flavor case by making the mass of {y^ — u T ) pair much smaller than 1 eV. 



14 



(This case will correspond to an inverted mass hierarchy p6|.) 

3 Cosmological Implications 

In studies of large scale structure, it has become clear that a critical den- 
sity Cold Dark Matter (CDM) universe with density perturbations that have 
roughly a Harrison- Zeldovich spectrum cannot simultaneous fit observations 
of structure on large and small scales. Normalized on large scales to fit the 
COBE observations, CDM produces too many clusters of galaxies and galac- 
tic pairwise velocities far in excess of observations ||47| . This basic problem 



is neatly resolved in a cold plus hot dark matter scenario because the growth 
of the density perturbations on small scales is damped by the presence of the 
hot dark matter component [7], p| . Since the neutrino oscillation experi- 
ments may be telling us that the neutrino component is spread among several 
flavors we would like to understand, in a quantitative way, how this changes 
the usual (one flavor) cold plus hot dark matter model. We will begin with a 
general description of the effects of a hot component on the growth of density 
fluctuations. First we will give some analytic formulae for estimating the size 
of the effects, then we will show some results of more detailed calculations. 
Finally we conclude with some remarks about the uncertainties in cosmolog- 
ical parameters which currently prevent one to use cosmological data to pin 
down the values of the neutrino masses. 
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3.1 Free— Streaming or Jeans Masses 

Neutrinos retain large velocities from the time before primordial nucleosyn- 
thesis when they were in thermal equilibrium. This means that if one initially 
sets up density perturbations, the neutrinos will rearrange themselves to a 
different pattern at some later time, erasing neutrino density fluctuations 
on length scales over which the neutrinos can have traveled. This "free- 
streaming" or "Jeans" length scale Aj is given at time t by 

Xj(t) = a(t)r(t)v(t), (6) 

where a(t) is the scale factor which describes the expansion of the universe, 
r(t) is the conformal time (dr = dt/a(t)), so a(t)r(t) is the physical horizon 
size, and v(t) is the average neutrino velocity. An analytic fit for v is given 
by PJ 



v(t) 



a - 2 



l+( — 

(X nr 



1-1/2 



(7) 



where a nr is the scale factor when the neutrinos become non-relativistic, and 
we use units where the velocity of light is unity. We have 

a nr = 2.25T 7 °/m„ (8) 

where T° = 2.35 x 10~ 4 eV is the present cosmic background photon tem- 
perature. The value of the Jeans' length at the current time to is 

Aj(t ) = St a nr = 3.2 ( — ) h- 1 Mpc. (9) 

h is the parametrization of our ignorance of the true value of the Hubble 
constant H = 100 h km s _1 Mpc -1 . To relate this length to the masses of 
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known astrophysical structures we convert to the total mass contained in the 
volume 47r(Aj/2) 3 /3, (the Jeans' mass). Currently, the Jeans mass Mj v (to) 
is 

Mj„(t ) = M Ho (a nr ) 3 = 4.6 x 10 12 (— V h~ l M , (10) 

where M Ho = 3.14 x 10 22 h' 1 M is the mass contained in the current Hubble 
volume. Interestingly enough, the mass in (10) is comparable to a galactic 
mass for neutrinos in the range m v ~ 1 — 4 eV. Note that in the formulae 
presented here we assume a critical density of matter, and only a fraction 
is in the form of hot dark matter- the remainder is comprised of cold dark 
matter and baryons. 

During the matter dominated epoch, p oc a~ 3 , v oc a" 1 (while the neu- 
trinos are non-relativistic) and ar oc a 3 / 2 , so the Jeans' mass Mj u oc a~ 3 / 2 
and decreases with time. When the neutrinos are relativistic, v ~ 1, and the 
Jeans' mass increases with time. The maximum in the Jeans mass occurs 
when the neutrinos are just becoming non-relativistic, i.e., when a(t) = a nr : 

/leV\ 3/2 

M^max) ~ M Ho (a nr f' 2 = 4 x 10 17 h' 1 M . (11) 

\ m v J 

This relation is approximate, because the time when neutrinos become rel- 
ativistic is quite close to the time when matter begins to dominate over 
radiation in the universe, so the horizon size does not scale exactly as a~ 3 / 2 . 
The exact formula must be determined numerically for each value of the neu- 
trino mass. The formula in (11) is accurate to within an order of magnitude, 
which is good enough for estimating the size of the effects under study here. 
In the older models of structure formation in which the dark matter is to- 
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tally composed of neutrinos of mass ~ 30 eV, this maximum Jeans mass is 
somewhat larger than a typical mass of a cluster of galaxies. 

In studies of single massive flavor C+HDM models, a more typical value 
of the neutrino mass ~ 6 eV (for h = 0.5) was preferred, yielding a maximum 
Jeans' mass of Mj v (max) = 3 x 10 16 h^ 1 M , which is slightly smaller than 
the so-called "great attractor" ES| , the largest concentration of mass which is 
pulling galaxies gravitationally on scales of order 40 h~ x Mpc. If we assume, 
as in scenario i), that the same Q u is now composed of 3 nearly degenerate 
neutrinos of mass ~ 2 eV, Mj v {max) = 1 x 10 17 h^ 1 M & , which is somewhat 
larger than the "great attractor" mass. 

We will proceed to examine how these two mass scales (M Ju (t ) and 
Mj v {max)) are incorporated in predictions of structure formation. 

3.2 Growth of Density Perturbations 

In the early universe, when matter begins to dominate the energy density 
of the universe, the CDM and C+HDM mass fluctuation spectra are iden- 
tical. The only feature imprinted on the initial spectrum at this time is a 
decrease of amplitude on length scales smaller than the horizon at matter 
domination. These scales experience a smaller amount of growth during the 
radiation dominated era due to the effects of radiation pressure. After matter 
domination the fluctuation amplitudes on these scales can grow, driven by 
their own self-gravitation. If a volume of space contains a mass M which is 
slightly different than the average M for that size, and the dark matter is 
strictly CDM, the fluctuation amplitude will grow proportional to the scale 



18 



factor a 

5M M-M 

(M) = oca (12) 

Growth in the C+HDM model is more interesting. When a mass fluctu- 
ation contains a mass M < Mj v (t) at time t, the neutrinos will stream out 
of the fluctuation and change the local gravitational potential. This smaller 
potential also reduces the growth rate of the mass fluctuations to |], |50 



<) M (M) (xaP- p = \ ( -1 + J25 - 24fi v ) (13) 



M 4 

The fact that the Jeans mass decreases with time means that if M < 
Mj v at some early time, sometime in the future M > Mj v will be true. 
Once M > Mj v) the growth rate returns to the CDM growth rate, as the 
neutrinos will then effectively be "cold" on this mass scale. This leads to 
scale dependent changes in the mass fluctuation spectrum. The equations 
for the Jeans mass (pT|) and (fL0|), and the growth rate (0) also illustrate 
the effect of having multiple flavors of massive neutrinos with comparable 
masses. The Jeans mass depends on the individual neutrino mass, while the 
growth rate depends only on the sum of the neutrino masses, as 

n„ = h- 2 y (14) 

4^94eV 1 1 

So increasing the number of degenerate flavors means that the Jeans mass 
will be increased for a given value of Q u . 

The growth of fluctuations on mass scale M in the intermediate region 
Mj v {max) < M < Mj v (now) will be "damped" relative to their value in 
a universe dominated only by CDM. The amount of damping between the 
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time when the neutrinos become non-relativistic and when the Jeans mass 

becomes smaller than scale M is given by 

Mj v (max) 



damping factor 



(2p-2)/3 

(15) 



M 

The maximum amount of damping occurs for mass scales M -C Mj v (to) 
which is a constant given by 

Mj v [max) 



max damping factor 



(2p-2)/3 

= (a nr ) 1 - p (16) 



M Jv (t ) 

The difference, then, between having the mass equally spread among M v 
flavors as opposed to 1 flavor results in an additional small scale damping of 
7VJ- p , which is between 10% - 20% for 0.20 < Vt p < 0.30. 

The combination of hot and cold components of dark matter thus has the 
following effects, when compared to the a model with only a cold component. 
First of all, the growth of density perturbations with masses in the interme- 
diate range Mj(now) < M < Mj(max) is reduced by the factor given in 



eq. fllq) which depends mainly on Vt v . Below Mj(now) the mass fluctua- 
tion spectrum follows the scale dependence of a CDM model but with an 
amplitude reduced by the factor in eq. fll6|) . Splitting the neutrino density 
among H v flavors of approximately degenerate neutrinos produces an addi- 
tional damping of N„~ l below a somewhat higher value of Mj u (max). This 
additional damping may be important for reconciling observations of clusters 
of galaxies with a COBE normalized spectrum of density fluctuations. 

3.3 Detailed Comparisons to Structure Formation 

In order to do a more careful comparison of these models to observed large 
scale structure, a more accurate evaluation of the effects expressed above 
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is required. We achieve this by integrating the linearized general relativis- 
tic evolution equations for the photons, neutrinos (massive and effectively 
massless flavors), baryons and CDM particles. The procedure we use is 
described in detail in ref. |5"If . The models have a baryon mass fraction 
Qbaryon = 0.0125//i 2 , which is reasonably consistent with nucleosynthesis 
[p2|. The CDM fraction is then given by Qcdm — 1 ~~ &v ~ ^baryon (since 
here we are only considering critical density universes). In figures 1 and 2, 
we use a Harrison-Zeldovich (scale free) spectrum normalized to COBE J53 



with h = 0.5. After smoothing the density field with a low pass Gaussian 
filter, we arrive at the curves shown in figures 1 and 2. The masses on the 
x-axis are (2ii) 3 ^ 2 r^p, as appropriate for a Gaussian filter radius of rf. In 
the three panels a), b), and c), we show curves for Q u = 0.20, 0.25, and 0.30 
respectively. Within each panel we show curves for 1, 2, and 3, degenerate 
neutrino flavors. In panel a), we also display the analogous CDM model 
curve for comparison. 

The models under consideration are now known to fit a large amount of 
data on large scale structure || [T2|, |TB|, ^0], but here we will concentrate 
on the two toughest constraints for C+HDM models to satisfy. The first 
constraint is on the amplitude of mass fluctuations which form clusters. It 
has been argued |54| that the number density of rich clusters requires that on 
a scale ~ 1O 15 M , 5M/M = 0.58 ±0.1 (2 times their "error bar"). As stated 



in ref. [54] the errors are hard to estimate because of possible systematic 
uncertainties in the masses assigned to clusters. Indeed, the mass scales for 



two different determinations of the cluster mass function in refs [55 and 
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[ 56[j differ by almost a factor of 2. Ref. |54| argues that the cluster masses 
could be overestimated, while measurements of the gravitational lensing by 
the clusters suggest [|59] that the usual mass estimates (used to determine 
(8M/M)) may be low by as much as a factor of 3. In order to allow for this 
uncertainty, we also give an estimated horizontal error bar which is a factor 
of two (in cluster mass) from the value 1 x 10 15 h' 1 M , to reflect the size 
of the possible systematic error. 

At the far left end of the figures we have also drawn a line which corre- 
sponds to the lower limit of the mass fluctuation amplitude implied by the 
measurements of the abundance of "damped Lyman alpha systems" as indi- 



cated by the recent survey of ref. ]60| ; see also ref. . While early work |J2 



indicated that C+HDM models (with Q u = 0.3) predictions were far below 
the observations, it now seems as though the assumptions used may have been 
too restrictive. There is little information about the exact nature of these 
damped Lyman alpha systems, other than their neutral hydrogen column 
densities and spectral line widths. A more general analysis |53] shows how 
to properly estimate number abundances of objects on these scales p3. We 
have used this latter procedure to estimate the amplitude of the constraint 
line in figures 1 and 2. It is still uncertain what mass is to be identified with 
these damped Lyman alpha systems. It seems reasonable that we should 
associate dark matter halo masses of between ~ 10 10 — 10 11 h^ 1 M Q and we 
have placed our constraint line there in our figures. (The curves should pass 
above the constraint line somewhere in the range.) The mixed dark matter 
models are more compatible with the lower end of this mass range. However, 



22 



with Vt v = 0.30 and h = 0.5 one has to take very small halo masses even 
to make this model work marginally. From figure 1 then it is apparent that 
Q u < 0.30 if h = 0.5 with a pure Harrison-Zeldovich spectrum. 

Figure 1 also points out an interesting aspect of spreading the neutrino 
mass density among a number of massive neutrino flavors H v . With M v = 1, 
the amplitude of the density fluctuations on cluster scales seems to be slightly 
larger than what is required; increasing M v improves the fit. This result 
depends only weakly on Q u on this mass scale. [Eq . fllBl) predicts less than 10% 
difference in the cluster scale 5M/M amplitude over the range Q u ~ 0.2—0.3.] 
However, because of the Q u dependence in the Ml~ x damping factor (viap), 
M v = 2 works best for Q u = 0.30 and M v = 3 works best for Q u = 0.20. 
On the other hand increasing the number of flavors makes the disagreement 
with Lyman alpha systems datum worse for Q u = 0.30, so M v > 1 seems to 
work best for models with Q u < 0.30. 

We should point out here that the conclusions we draw from figure 1 
hold only for the spectrum, Hubble constant and baryon fraction used in the 
models, and altering any of these parameters changes the conclusions. We 
will explore this further later. The purpose of figure 1 is mainly to illustrate 
the effects of modifying the hot dark matter composition. 

3.4 Additional Models Which Can Accommodate the 
LSND Results 

In the above we have considered the consequences of the three neutrino sce- 
narios in which we have H v flavors with roughly equal mass. However, if we 
take the LSND results as a serious constraint, then scenario (iii) admits the 
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possibility that all three known neutrino flavors can have masses in the eV 
range, although the electron neutrino mass can still be significantly different 
from the fj, and r neutrino masses (which form a degenerate pair). Since only 
the square of the mass difference is measured, we have two cases to consider: 
a) m Vtl > m Ue , ("2+1") and b) m Ue > m V)1 ("1+2"). 

In a "2 + 1" model, where and v T are each heavier than v e , the lightest 
neutrino (v e ) will have a very large maximum Jeans' mass, (M l Jiy (max)) , com- 
pared to the other two heavier neutrino Jeans masses (Mj v (max)). In the 
"1 + 2" model, wherez/^ and v T are lighter than v e , the lightest neutrinos (v^ 
and v T ) will have very large M l Jv (max), compared to the remaining heavier 
v e neutrino Mj (max). For M l Jv (max) > M > Mj v (max), there will be an 
amount of damping which would be there if only the lightest neutrino com- 
prised VL U . For Mj v (max) > M the damping will be intermediate between 
the one and three flavor massive neutrino case. 

To illustrate these two cases, we will use values of Am 2 = 8 and 20 eV 2 , 
for the "2 + 1" and "1 + 2" models respectively in the context of an Q u = 0.25, 
h = 0.5, n — 1 cosmology. In figure 2 we plot the curves for the 2 + 1 model 
(m Ue = 0.22 eV and m Vr ,m Vii = 2.84 eV), and the 1 + 2 model (m Ve = 4.524 
eV and m Ur , m v = 0.685 eV. For comparison we also plot the curves for 
M v = 1, 2, and 3 flavors for the same Vl u (= 0.25). The curve for 2 + 1 (1 + 2) 
lies in between the curves for the 2 and 3 (1 and 2) degenerate flavors. 

We can estimate the difference between having 3 flavors with equal and 
unequal masses as follows. In the "1+2" and "2+1" models the lightest 
flavor mass m v \ is lighter than the three equal flavor mass m u3 . The heavier 
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flavor mass m v h is greater than m v %. Consider a mass scale M < Mj (max) . 

After some time, when Mj (t) = M, the amount of damping is the same 

in both the equal and unequal mass cases. After that time that scale will 

experience an additional damping of \M l Jv (max) / 'Mj v (max)]( 2pi ~ 2 ^ 3 (eq. 

in the unequal mass case, and [Mj 3 ^ lavors \max) /Mj^(max)]( 2p ~ 2 ^ 3 in the 

3 flavor case. Here pi is p in equation ( O ) calculated assuming Q u is only 

made from the lightest neutrino(s). The ratio of the damping in the 1+2 

(and 2+1) cases to the 3 equal mass case is (for scales < Mj v (max)) 

, . .. (mui/mvh) {1 - pi) n7 , 
damping ratio = ; -7- — r, (17) 

(m^/m^-py 

The ratio in ( \L7\ j is always larger than 1, indicating that the damping is 
greatest in the 3 equal flavor case. However, if the ratio of the masses of 
heavy to light flavors is < 4 the difference between damping in the various 
cases differs from the 3 flavor case by < 2 %. This is the reason for choosing 
the particular Am 2 for plotting. Had we plotted, for example, the 2+1 and 
1+2 models using Am 2 = 6 eV 2 , (the same value used in ref. |4j], but for 
Q u = 0.20) both of the curves would be nearly indistinguishable from the 
3 degenerate neutrino case for masses < Mj (max). When the light and 
heavy masses differ by an order of magnitude or more, as in figure 2, the 
2 + 1 model goes to the 2 degenerate flavor model and the 1 + 2 goes to 
the single flavor model. We note that because the lightest neutrino weighs 
less than the 3 degenerate flavor neutrino, on the very largest scales where 
M l Jv (max) < M < Mj v ^ lavors \max) , there will be some damping which is 
absent in the 3 degenerate flavor case. 

We see that these additional models offer now a continuum of AM/M 
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values between the 1 and 3 degenerate mass neutrino cases. The sole mo- 
tivation for considering where the neutrinos have masses and mass 
differences of order a few eV in scenario iii) is the possibility of detecting 
oscillations with Am 2 in the few eV 2 range in the LSND (and also KAR- 
MEN) experiments. In the framework of GUT theories however, such a mass 
spectrum will not be easy to understand. 

3.5 Sensitivity to Cosmological Parameters 

In the previous sections we have plotted results for specific values of the 
Hubble constant, initial density fluctuation spectrum, and flbaryon- We have 
not drawn strong conclusions about the neutrino mass scenarios, because the 
results are quite sensitive to the specific values used for these cosmological 
parameters. We will discuss the effects of these parameters on the cosmolog- 
ical structure formation and give specific examples of models from each of 
the three mass scenarios. 

1. Hq: The value of the Hubble constant has been the subject of a long 



and ongoing observational campaign (see [65] for some recent measure- 
ments). The experiments find that 0.4 < h < 1.0. This uncertainty is 
quite large when seen from the point of view of the required neutrino 



masses. Rewriting equation (|14| ) for degenerate neutrinos we find 

-= 4 4^fe) 2eV - (18) 

which emphasizes that the neutrino masses scale as h 2 . In figures 1 and 
2 we have used h = 0.5. Some of the recent measurements [pql seem 
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to imply a larger value of h. In particular, if h = 0.6, then the masses 
used in the previous examples are increased by a factor of 1.44. For 
example, ref. M advocates using Af u = 2 and Q u = 0.2 with Am 2 = 6 
eV 2 , h = 0.5. If instead h = 0.6 with the other parameters fixed, then 
the oscillation signal at Los Alamos would be consistent with Am 2 = 12 
eV 2 . 

The scale which corresponds to the horizon size at matter domination 
also depends on the Hubble constant. If we use a larger value of h, then 
matter domination occurs earlier, when the horizon size was smaller, 
which means that the radiation dominated era is less effective at arrest- 
ing the growth of small scale fluctuations. Such an effect exacerbates 
the small scale problems in CDM models and favors using a larger value 

of a, §,0. 

are attributed 

2. Initial Mass Fluctuation Spectrum 

It has long been known that inflation predicts a power spectrum of 
density fluctuations with a spectral index n close to unity. Since the 
amount of deviation is strongly model dependent, many investigators 
are content to use n — 1. In figures 1 and 2 we have also used the value 
of n — 1 for the spectral index. However, even small deviations from 
n = 1 can lead to significant changes in the conclusions. To see this we 
give the initial spectrum in terms of rms mass fluctuations (^f)' 

AM 
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M -(3+n)/6 ( 1Q ) 



In the simplest models of inflation, particularly those based on GUTs 
(see, e.g., ref. [|9|, |30|), n ~ 0.94 — 0.98, although other values of 
n are certainly possible |66[. From eq. ( |T9"D we see that decreasing 



(increasing) n from unity decreases (increases) the small scale power. 
Since the models are normalized at large scales to COBE observations 
the mass fluctuation curves in figures 1 and 2 will "pivot" around a 
very large mass scale ~ 10 21 h~ x M . The "pivot" mass scale here is 
somewhat smaller than the horizon mass, because the best fit COBE 
quadrupole anisotropy scales as e 1_n f53fl . (We have fit our spectra by 
normalizing to the 7th multipole moment of the Sachs- Wolfe anisotropy, 
as recommended in ref. [|53] ) 



3. Qbaryon- In the past few years comparisons of the primordial light ele- 
ment abundances inferred from observations with those calculated have 
led to strict limits |52|] on the amount of baryonic material in the uni- 
verse, 0.04 < Q baryon (0.5/h) 2 < 0.06. We have used the central value 
of this range for the models in figures 1 and 2. Improvements in the 
observations of deuterium and 3 He and determinations of the neutron 
lifetime have now led to a situation where all of the light element abun- 
dances are not consistent in a universe where the three known neutrino 
flavors have the standard number densities (derived from thermal equi- 



librium) at nucleosynthesis |p7 |. If, however, the He abundances de- 



rived from observations have been systematically underestimated |)8 
by about 5%-10% then big bang nucleosynthesis would be made con- 
sistent with the three flavors of neutrinos, provided that the baryon 
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density is about 50% larger than previously thought. This would go in 
the right direction to explain why the baryon to dark matter ratio in 
cluster cores is so much larger than flbaryon [0 ■ 

If we allow for larger baryon fractions, this will also change the am- 
plitude of the mass fluctuations for a given Q u and n. We can un- 
derstand this through the following. Baryons (mostly protons) after 
nucleosynthesis are electrically charged and so are strongly coupled to 
the photon field. Baryonic density fluctuations cannot grow until the 
photon temperature decreases to allow the stable formation of neutral 
hydrogen. Because photon-baryon decoupling does not happen until 
after matter domination, the baryons are prevented from falling into 
the gravitational wells supplied by the cold dark matter. This leads to 
a damping in the growth of density fluctuations (relative to with 
no baryons). The amount of damping (on scales smaller than the hori- 
zon size at photon baryon decoupling) in the final amplitude of density 
fluctuations is constant, so the effect is quite similar to changing the 
number of degenerate neutrino flavors. 

These uncertainties mean that one cannot currently use the cosmological 
data to determine which of our three scenarios is correct, or even what the 
value of VL V is. One must do a systematic study of the available parameter 
space taking into account the full range of these parameters. Studies of the 
interplay between these cosmological parameters have been done in various 



contexts [13. 29 , 7C], with more in progress [71 . 
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As a graphic illustration of this uncertainty, and also to demonstrate use 
of the formulae presented in sections 3.1-3.3, we show four models in figure 
3 that give nearly identical predictions of structure formation, but use quite 
different solutions of the neutrino oscillation data. The choices of parameters 
are as follows. 

1. The solid curve. [Scenario (iii)]. This model with M v = 2, Q u = 0.20, 
n = l, and h = 0.5 is advocated as a good fit to cosmological obser- 
vations by ref. f|] consistent with the Los Alamos LSND experiment. 
We present this model for comparison. 

2. The dot-dashed curve. [Also scenario (iii)]. We increase h to 0.6 in 
this M v = 2 model, but now there will be too much small scale power. 
To compensate for this, we increase Q u to 0.25, and decrease n to 
0.95, a value which actually is more in line with standard inflationary 
predictions. The v e — experiments would then be expected to see a 
signal corresponding to Am 2 = 18 eV 2 . The baryon fraction has been 
scaled from the first model as h 2 consistent with the nucleosynthesis 
constraints. 

3. The short dashed curve. [Scenario (i)]. Again using h = 0.6 and 

we now consider M u = 3 degenerate flavors. Because the extra flavor 
provides an additional damping factor over the two flavor model of 
(2/3) 16 = 0.94, (eq. [16]), to get a curve similar to the previous one 
we can increase n to 0.98, which increases the amplitude at the cluster 
mass by (10 21 /10 15 ) ' 0075 = 1.10, (eq. [19|), so we get very nearly the 



same fit at the cluster scale. The value n = 0.98 happens to be the 
value predicted in a particular SUSY GUT inflation model pOfl . 

4. The long dashed curve. [Scenario (ii)]. Here we return to h = 0.5, 
Vt v = 0.2, and n = 1 as in the solid curve. To compensate for the 
increase in amplitude by going back to 1 flavor (2/l) a13 = 1.09 in 
the solid curve, we increase flbaryon to 0.10. This value of Qbaryon is 
consistent with 3 relativistic neutrino flavors during nucleosynthesis if 
the 4 He abundances have been systematically underestimated. 

We have now shown that allowing for more than one degenerate (in 
mass) neutrino state, as indicated by solar and atmospheric oscillation exper- 
iments, produces a "degeneracy" in the predictions of structure formation for 
C+HDM models, given the uncertainty in cosmological parameters. In order 
to break this latter "degeneracy" we need a convincing detection of neutrino 
oscillations in an accelerator experiment, or an improvement in the determi- 
nation of cosmological parameters, or both. In the meantime, we just note 
that there is a very rich structure contained in C+HDM models of structure 
formation, which still allows considerable flexibility in fitting astrophysical 
data. 

4 Conclusions 

We have investigated the impact of neutrino oscillations, as indicated by a 
number of solar and atmospheric neutrino experiments, on the 'cold plus 
hot' dark matter scenario of large scale structure formation. We are led 
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to three distinct scenarios for neutrino masses and mixings with interesting 
predictions for the ongoing/planned experiments. We note in particular the 
expectations for — u e and — u T oscillations being currently searched for. 
The cosmological implications of the three scenarios are explored in some 
detail. For some choices of the cosmological parameters (particularly h = 0.5, 
n — 1, and Vt haryon = 0.05), the two and three neutrino flavors scenarios 
provide a somewhat better fit to the present data than the single (neutrino) 
flavor case. However, as we show, this need not hold for a different parameter 
choice. It is too early to pick out the best model but it is remarkable that 
taking account of the oscillation experiments, the C+HDM models can still 
provide a good fit to the large scale structure data within the context of the 
simplest inflation models. 
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5 Figure Captions 

Fig. 1. We show the difference in the rms mass fluctuation amplitude 
(AM/M) between having the hot dark matter (neutrino) component dis- 
tributed among 1, 2, or 3 flavors of degenerate neutrinos which correspond 
to our scenarios ii), iii) and i) respectively, for fitting the neutrino oscillation 
data. Here h = 0.5, Qbaryon = 0.05, and n — 1, with the spectrum normal- 
ized to COBE. The data points are: 1) the mass fluctuation estimate based 



on cluster formation |54 and 2) the lower limit on the mass fluctuations 
consistent with the formation of Lyman alpha cloud abundances at high red- 
shift (The mass at which this lower limit is to be imposed should be 
somewhere along the constraint line.) Panels a), b), and c) correspond to 
Q u = 0.20, 0.25 and 0.30, respectively. It can be seen that increasing the 
number of neutrino flavors marginally improves the fit to the cluster data. 
In panel a) we also show the curve for a pure cold dark matter model. 

Fig. 2. The rms mass fluctuation amplitude in a special case of neutrino 
scenario iii, which has been made consistent with the LSND best fit of 5m? = 
6 eV 2 in a model with Q u = 0.25, h = 0.5, n — 1. For comparison we also 
show (from top to bottom) the same f2„ with 1, 2, or 3 degenerate neutrino 
flavors. The "2 + 1" and "1 + 2" model curves are intermediate to the 
degenerate neutrino cases. The data are the same as in Figure 1. 

Fig. 3. The rms mass fluctuation in some selected models with varying 
cosmological parameters. We plot curves for two models from scenario iii 
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with M v = 2 for (fi„ = 0.20, h = 0.5, n = 1) - as advocated in ref. - 
and for = 0.25, h = 0.6, n = 0.95). These two models would predict 
Sm 2 = 5.5 and 18 eV 2 at the LSND and KARMEN experiments. The third 
curve is for a model with J\T V = 3 (Q v = 0.25, h = 0.6, n = 0.98). Lastly we 
show a curve from scenario ii) with M v (Q u = 0.20, h = 0.5, n = 1.00). Here 
we have increased O. haryon to reduce the ampitude at cluster scales. As can be 
seen, it is difficult to determine which neutrino mass scenario is correct based 
on current cosmological data alone. The data are the same as in Figure 1. 
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Figure 1 . Mass Fluctuation Amplitude for Degenerate Neutrinos 
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Figure 2. Mass Fluctuation Amplitude 




Figure 3. The Effect of Cosmological Parameter Uncertainties 
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